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Abstract 

Neutrino scattering on nucleons in the regime of deeply virtual kinematics is studied both 
in the charged and the neutral electroweak sectors using a formalism developed by Blumlein, 
Robaschik, Geyer and Collaborators for the analysis of the Virtual Compton amplitude in 
the generalized Bjorken region. We discuss the structure of the leading twist amplitudes of 
the process. 



Exclusive processes mediated by the weak force are an area of investigation which may gather a 
wide interest in the forthcoming years due to the various experimental proposals to detect neutrino 
oscillations at intermediate energy using neutrino factories and superbeams 0J. These proposals 
require a study of the neutrino-nucleon interaction over a wide range of energy starting from the 
elastic/quasi-elastic domain up to the deep inelastic scattering (DIS) region (see [2j,[3j,[4j,[5j,[6j 
for an overall overview). However, the discussion of the neutrino nucleon interaction has, so far, 
been confined either to the DIS region or to the form factor/nucleon resonance region, while the 
intermediate energy region, at this time, remains unexplored also theoretically. Clearly, to achieve 
a "continuos" description of the underlying strong interaction dynamics, from the resonant to the 
perturbative regime, will require considerable effort, since it is experimentally and theoretically 
difficult to disentangle a perturbative from a non-perturbative dynamics at intermediate energy, 
which appear to be superimposed. This is best exemplified - at least in the case of electromagnetic 
processes, such as Compton scattering - in the dependence of the intermediate energy description 
on the momentum transfer 0. In this respect, the interaction of neutrinos with the constituents 
of the nucleon is no different, once the partonic structure of the target is resolved. From our 
viewpoint, the presence of such a gap in our knowledge well justifies any attempt to improve the 
current situation. 

Together with Amore, we have pointed out |Hj that exclusive processes of DVCS-type (Deeply 
Virtual Compton Scattering) could be relevant also in the theoretical study of the exclusive neu- 
trino/nucleon interaction. Thanks to the presence of an on-shell photon emitted in the final state, 
this particle could be tagged together with the recoiling nucleon in a large underground detector 
in order to trigger on the process and exclude contamination from other backgrounds. With these 
motivations, a study of the vN — > vN^ process has been performed in [HJ. The process is mediated 
by a neutral current and is particularly clean since there is no Bethe-Heitler contribution. It has 
been termed Deeply Virtual Neutrino Scattering or DVNS and requires in its partonic description 
the electroweak analogue of the "non-forward parton distributions", previously introduced in the 
study of DVCS. 

In this work we extend that analysis and provide, in part, a generalization of those results to 
the charged current case. Our treatment, here, is purposely short. The method that we use for 
the study of the charged processes is based on the formalism of the non-local operator product 
expansion and the technique of the harmonic polynomials, which allows to classify the various 
contributions to the interaction in terms of operators of a definite geometrical twist |TT]. We 
present here a classification of the leading twist amplitudes of the charged process while a detailed 
phenomenological analysis useful for future experimental searches will be given elsewhere. 

2 The Generalized Bjorken Region and DVCS 

Fig. [U illustrates the process that we are going to study, where a neutrino of momentum I scatters 
off a nucleon of momentum P\ via a neutral or a charged current interaction; from the final state 
a photon and a nucleon emerge, of momenta q2 and P2 respectively, while the momenta of the 
final lepton is I'. We recall that Compton scattering has been investigated in the near past by 
several groups, since the original works [TJJ [THJ [TJj . A previous study of the Virtual Compton 
process in the generalized Bjorken region, of which DVCS is just a particular case, can be found 
in [in]. From the hadronic side, the description of the interaction proceeds via new constructs 
of the parton model termed generalized parton distributions (GPD) or also non-forward parton 
distributions. The kinematics for the study of GPD's is characterized by a deep virtuality of the 
exchanged photon in the initial interaction (u+p — > u+p + ^) ( Q 2 m 2 GeV 2 ), with the final state 
photon kept on-shell; large energy of the hadronic system {W 2 > 6 GeV 2 ) above the resonance 
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Figure 1: Leading hand-bag diagrams for the process 



domain and small momentum transfers \t\ < 1 GeV 2 . In the electroweak case, photon emission 
can occur from the final state electron (in the case of charged current interactions) and provides 
an additional contribution to the virtual Compton amplitude. We choose symmetric defining 
momenta and use as independent variables the average of the hadron and gauge bosons momenta 

Pl,2 = PTj ?1,2=?±|, (1) 

with A = P2 — Pi being the momentum transfer. Clearly 

A 2 

P ■ A = 0, P 2 = M 2 - — (2) 

and M is the nucleon mass. There are two scaling variables which are identified in the process, 
since 3 scalar products can grow large in the generalized Bjorken limit: g 2 , A • q, P ■ q. 

The momentum transfer A 2 is a small parameter in the process. Momentum asymmetries 
between the initial and the final state nucleon are measured by two scaling parameters, £ and rj, 
related to ratios of the former invariants 



^ = -2p-q ^=2p-q ® 

where £ is a variable of Bjorken type, expressed in terms of average momenta rather than 
nucleon and gauge bosons momenta. The standard Bjorken variable x = —q\j(2P\ ■ q\) is trivially 
related to £ in the t — limit and in the DVCS case r] = — £. 

Notice also that the parameter £ measures the ratio between the plus component of the mo- 
mentum transfer and the average momentum. 

£, therefore, parametrizes the large component of the momentum transfer A, which can be 
generically described as 

A = -2£P - A ± (4) 



where all the components of A^ are O (^y |A 2 

3 Bethe-Heitler Contributions 

Prior to embark on the discussion of the virtual Compton contribution, we quote the result for 
the Bethe-Heitler (BH) subprocess, which makes its first appearance in the charged current case, 
since a real photon can be radiated off the leg of the final state lepton. The amplitude of the BH 
contribution for a W + exchange is as follows 



T ^ = -\e\-^=-^u{l') 



[i — Ay + it 



2 



where e is the polarization vector of the photon and 
ii0 9 - 



T, 



w~ 

BH 



U{P 2 



v{l) 



(l - A) 2 + ie 



A 2 — + 



(i^A 2 ) - i^(A 2 )) 7 5 + (F 2 "(A 2 ) - F 2 d (A 2v 



a 



2M 



U{P 1 



(5) 



(6) 



for the W case, with D u& \qi) / (A 2 — + ie) being the propagator of the W's and F 12 the usual 
nucleon form factors (see also [8j). 



4 Structure of the Compton amplitude for charged and neu- 
tral currents 

Moving to the Compton amplitude for charged and neutral currents, this can be expressed in 
terms of the correlator of currents 



Tfxu i 



d 4 xe^(P 2 \T (jJ(;r/2) ± > z °{-x/2)) \Pi) , 
where for the charged and neutral currents we have the following expressions 



(7) 



J^ z °(-x/2) = 
r w+ (-x/2) 

r w ~{-x/2) 



9 T 



2 cos 9 W 



^ u {-x/2)^{g z uV + g^ uAl ^ u (-x/2) + ^ d (-x/2)^(g* v + ^ A7 & )^(-a;/2) 



„,5 



U ^(-x/2) 7 ^(l-7 5 )^(-x/2), 



2^2 
9 

2V2 



^ d (-x/2)r(l-^)U du M-x/2), 



J^(x/2) = ^ d (x/2)Y (-ie) + rl>«(x/2)-f (§e) ^(x/2) . 

Here we have chosen a simple representation of the flavour mixing matrix U* d = U u d = Udu 
cos 9c, where 6c is the Cabibbo angle. 
The coefficients gy and g\ are 
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9uv 



9dv 



+ - sin 2 8 
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9d = ~ 
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(9) 



(10) 



are the absolute values of the charges of the up and down quarks in units of the electron charge. 
A short computation gives 



{P 2 \T(.r v {x/2)J z \-x/2))\P l ) = 
(PS u (x/2)g ulu S(x) lfl (g z v + g z uAl ^ u {-x/2) 
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^ d {x/2) llx {g z dV + g z dA ^)S(-x)g dl ^ d ^/nPi), 

(11) 

{P 2 \T(j2(x/2)J^\-x/2))\P 1 ) = 
(P 2 \^(-x/2) lfi (l - ^)U ud S{-x) lu {-g d ) i/> d (x/2) + 
^ u {x/2) lv (g u ) S(x)^(l - ^)U ud M-*m\Pi), 

(12) 

(P 2 \T(rj(x/2)J^(-x/2))\P 1 } = 

(P 2 \ - ^ d {x/2)g dlv S{x) llx {\ - 1 b )U du ^ u {-x/2) + 

^ d (-x/2)^(l - 7 5 )5(-x)f/ du ^^/2)|Pi) , (13) 

where all the factors g/2\/2 and g/2cos9w, for semplicity, have been suppressed and we have 
defined 



S u (x) = S d (x) 



2tt 2 (x 2 — ie) 2 



(14) 



Using the following identities 



5 _ c P 5 _ • (i 



(15) 



we rewrite the correlators as 



Tf? = i / d A x 



- /IV 



2tt 2 (x 2 — ie) 



\9u9uV [ 



~9d9dv {S^ctvpO^ — ie^avpO^f^J + gd9dA 

(16) 



- /J.U 



rpW 
[IV 



T 

27r 2 (x 2 — ie) 2 



d 4 x 



iS, 
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\Pl), 



iqx ajj 
■ i ;4 ° ^ du 

t ax- 



2n 2 (x 2 - ief 



~iS, av p (<5£ + O b £) - e, avP (Oi + Of)] \P X 



(17) 



(18) 



We have suppressed the x-dependence of the operators in the former equations. The relevant 
operators are denoted by 

Of (x/2, -x/2) = ^(x/2) 7 ^ (-x/2) + ^ a (-x/2)^ a (x/2), 
6f (x/2, -x/2) =j0 a (x/2) 7 5 7 ^a(-^/2) r f a (-x/2) 7 Vi(x/2), 
Of (x/2, -x/2) = ^(x/2) 7 > a (-x/2) - i/> a (-x/2)<fi/, a (x/2), 
Of '(x/2, -x/2) = V> a (*/2) 7 5 yVa(-z/2) + t(-x/2) 7 5 7^(x/2), 



Of d (x/2, -x/2) = g u ^ u (x/2)^ d (-x/2) + g^ u {-x/2)^ d (^l 
0%(x/2, -x/2) = ^ tt (x/2) 7 V^(-^/2) - ? A(^/2)7 5 7%(x/2), 



(19) 
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(20) 



and similar ones with u <-> d interchanged. 

We use isospin symmetry to relate flavour nondiagonal operators {Off) to flavour diagonal 
ones {Off) 



(p\O ud {x)\n) = (p\O ud {x)r-\n) = (p| O ud (x),r- \n) = (p\O uu (x)\p) - {p\O dd {x)\p) 



udi 



\ddi 



{p\O ud (x)\n) = (n\0 M {x)\n) - (n\O uu {x)\n) 
(n\O du {x)\p) = {p\6 uu {x)\p) - (p\O dd {x)\p) , 



{n\O du {x)\p) = (n\O dd {x)\n) - (n\O uu {x)\n), 

where 

r = r ± r w 

are isospin raising/lowering operators expressed in terms of Pauli matrices. 



(21) 
(22) 



5 Parameterization of nonforward matrix elements 

The extraction of the leading twist contribution to the handbag diagram is performed using the 
geometrical twist expansion, as developed in [TOl 112 EH] , adapted to our case. We set the twist-2 
expansions on the light cone (with x 2 = 0) and we choose the light-cone gauge to remove the 
gauge link 



(P 2 \^ a {-kx) 1 ^a{kx)\Pl) 



tw.2 



J Dze-*t*- p *)F<* v \zi, z 2 , P { ■ P jX 2 , Pi ■ P j )U{P 2 ) [Y - ikP/tf ] U{P X ) + 
J Dze- ik ^G a ^{z h z 2 , P t ■ P jX 2 , P t ■ P j )U{P 2 



{ia» a A a ) 



ia ap x a A p 
M 



with < k < 1 a scalar parameter, with 

P z = P lZl + P 2 z 2 , 

and 



U{Px) , 
(23) 

(24) 



(p^-kxrfYMkxmY™- 2 = 

J Dze- ikix - p ^F 5aiv) {z u z 2 , Pi ■ Pjx 2 , Pi ■ Pj)U{P 2 ) [ 7 V - ikP/f$ ] U{P X ) + 
J Dze- ik ^G 5a ^{ Zl , z 2 , Pi ■ P.x 2 , Pi ■ P j )U{P 2 ) 1 h 



ikP 7 



ia a ^x a A p 
M 



U{Pi 



(25) 



The index {v) in the expressions of the distribution functions F, G has been introduced in order 
to distinguish them from the parameterization given in |10[ll6j. which are related to linear combi- 
nations of electromagnetic correlators. In the expressions above a is a flavour index and we have 
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and P2, and the measure of integration is defined by [lOj 



Dz = -d Zl dz 2 9(1 - zi) 9(1 + Zl ) 9(1 - z 2 ) 9(1 + z 2 ). (26) 

In our parameterization of the correlators we are omitting the so called "trace-terms" (see 
ref. since these terms vanish on shell. In order to arrive at a partonic interpretation one 

introduces variables z + and Z- conjugated to 2P and A and defined as 

Z + = 1/2(3! + 23), 

Z- = 1/2(^2 - Zl), 

Dz = dz + dz_9(l + z + + z_)0(l + z + - z_)0(l - z+ + z_)0(l - z+ - z_). (27) 

In terms of these new variables P 2 = 2Pz + + Az_, which will be used below. 

At this stage, we can proceed to calculate the hadronic tensor by performing the x-space 
integrations. This will be illustrated in the case of the W + current, the others being similar. We 
define 
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{P2\SfxavpO a I Pi) — S* v , 
(P2\e^auf30 a ^\P 1 ) = 



5a 



(P2\e,avp6 5a P\P l }=6 



~5a 



(28) 



and introduce the variables 



(29) 



where (z) is now meant to denote both variables (z + ,z_). The presence of a new variable Q 2 , 
compared to jTH|, is related to the fact that we are parameterizing each single bilinear covariant 
rather then linear combinations of them, as in the electromagnetic case. 
After some re-arrangements we get 



S ^ = ^lDz mr ^±{\ 



(Qi + »e) 



(Qj + ie 



F a( - U )(z) r 



+ Quip + %lv] + [Pz^lu + Pzvlp. 
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+ + %lv\ - [Pzulv + Pzvl^ 

Pz^Pzv + PziAv + Pzv% - 9nv{Pz ■ q) 



+ 
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M ^ M 

-PzuPzv + Pz^qu + Pzv% - g^u(P z ■ q) 



2 " M z 



M 



(30) 



with an analogous expressions for Sf£, that we omit, since it can be recovered by performing the 
substitutions 



7m- 7 5 7m ^-tV", 

pa(y) Qa(v) _^ p5a(u) q5ci(u) 



(31) 



in (TO . 

Similarly, for e" we get 



g d J DzF<»\z) 
g u J DzG«"\z) 
g d J DzG< v \z) 



+ ie) 
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{Ql + *e) 



-fiauP 



q a ^ - 



{Ql + ie) 



+ 



{Ql + ie 



~€[iavP 



(Ql + ie) W 



{Ql + ie) 
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M M{Ql + ie) 
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M M(Ql + ie) 



(32) 



The expression of can be obtained in a similar way. 



To compute the Tjf° tensor we need to include the following operators, which are related to 
the previous ones by g u , ga — > 1 



dx 



dx 



2n 2 (x 2 — ie) 2 
2tt 2 (x 2 — ie) 2 
2tt 2 (x 2 — ie) 



{P2\S IMXV pO aP \P 1 ) = 3- 
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{P2\S, av pO^\P 1 ) = 3 



5a 
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-{P 2 \e mv pO^\P x ) = ~el 



dx 



2tt 2 (x 2 — ie) 2 
In this CclSG el simple manipulation of (f30|) gives 



(P^pO^P,) 



fc (J,V 



(33) 
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{I- 



{Ql + ie 
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G a(»), 
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ia aP q a Ap 
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The expression of is obtained from (|34"|) by the replacements (|3TJ) . 
For the e a case, a re-arrangement of (f32j) gives 



(Qi + te 
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+ ie 
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{Ql + ie 
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Also in this case, the expression of the tensor is obtained by the replacements (f3*Tl) . 



(34) 



(35) 



6 The partonic interpretation 

At a first sight, the functions G^ u \ F 5( - u \ G 5 ^ do not have a simple partonic interpretation. 
To progress in this direction it is useful to perform the expansions of the propagators 



Qi + ie 2(P-q)[z+-£ + 7)Z- + ieY 

Ql + ie ^ ~2{P-q) [z+ + £ + riz- - ie] ^ ' 

which are valid only in the asymptotic limit. In this limit only the large kinematical invariants 
and their (finite) ratios are kept. In this expansion the physical scaling variable £ appears quite 
naturally and one is led to introduce a new linear combination 

t = z + + rjz- , (37) 

to obtain 

1 1 1 



Qi + ie 2{P . q )[t-Z + ie]' 

111 
Qf+Te ~ ~2(P-q)[t + £-ieY (38) 

Analogously, we rewrite P z using the variables {t, z_} 

P z = 2Pt + nz_ , (39) 
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7T = A + 2£P . 



(40) 



This quantity is strictly related to Aj_, as given in (J3J. The dominant (large) components of 
the process are related to the collinear contributions, and in our calculation the contributions 
proportional to the vector n will be neglected. This, of course, introduces a violation of the 
transversality of the process of 0(A±/2P ■ q). 

Adopting the new variables {t, zJ\ and the conjugate ones |P, 7r|, the relevant integrals that 
we need to "reduce" to a single (partonic) variable are contained in the expressions 
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(41) 



Here H(z + ,z_) is a generic symbol for any of the functions. We have similar expressions for the 
integrals depending on the momenta Q2. 

The integration over the z_ variable in the integrals shown above is performed by introducing 
a suitable spectral representation of the function H(t, +£z~, zJ). As shown in [10], we can classify 
these representations by the n = 0, 1, powers of the variable z_ , 



dz-z_ n h(- + f z_ 



(42) 



With the help of this relation one obtains 



t< 



dz.z„ n H(t + £z_,z_) 



1 dr 



r 



r n h n (t/r,0 



(43) 



The result of this manipulation is to generate single- valued distribution amplitudes from double- 
valued ones. In the single-valued distributions h n (t,^) the new scaling variables t and £ have a 
partonic interpretation. £ measures the asymmetry between the momenta of the initial and final 
states, while it can be checked that the support of the variable t is the interval [—1,1]. The twist-2 
part of the Compton amplitude is related only to the n = moment of Z-. Before performing the 
z_ integration in each integral of Eq. (jH| using Eq. ijlHJl - a typical example is Pq^(£) - we reduce 
such integrals to the sum of two terms using the identity 
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As shown in [TH], after the z_ integration, the integrals in (J4T1) can be re- written in the form 
1 



2P-q 



-1 (t - £ + ie) 
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(t-Z + ie) 



{{2q^P u + 2q v P" - g^2q ■ P)} + O^ir"), 



(45) 



where, again, we are neglecting contributions from the terms proportional to 7r M , subleading in 
the deeply virtual limit. The quantities that actually have a strict partonic interpretation are the 
^o(^0 functions, as argued in ref. [11]. The identification of the leading twist contributions is 
performed exactly as in |TQ1 - We use a suitable form of the polarization vectors (for the gauge 
bosons) to generate the helicity components of the amplitudes and perform the asymptotic (DVCS) 



limit in order to identify the leading terms. Terms of 0{\/ \J2P ■ q) are suppressed and are not 
kept into account. Below we will show only the tensor structures which survive after this limit. 



7 Organizing the Compton amplitudes 

In order to give a more compact expression for the amplitudes of our processes we define 

Tjxv 
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fj, pu u p/i 

9 + (q-P) + (q-P) 
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(46) 



Calculating all the integrals in the Eqs. (|30"|). (f3l|) and (|32|). (|35|) . we rewrite the expressions of 
the amplitudes as follows 
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where, suppressing all the subleading terms in the tensor structures, we get 
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M 



(48) 



while for the e afiu expression we obtain 



U(P 2 )e a " u U(P 1 



(2P ■ qf 



dt(3{t) [U{P 2 )i U(P 1 )fS(t,Z) + 

u{P2)ii a -^ l )u{p 1 )m,z) 



(49) 
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U{P 2 )S a ^U(P 1 



dt 



+ 



-i 



9 



2P ■ q 

l / 



V(p 2 )h u{p 1 )f^o + u{p 2 ){i 



M 



+ 



dt 



1 



-i 



IE 



t + i-ie {P-qf 



u(p 2 )i umtfsbt) + C/(P 2 )(i^|^)C/(Pi)^o (^0 



(50) 



U{P 2 )e a ^U{P 1 ) = f 1 dt( 1 1 ) e^ 2q ^\ x 

v ' K ' J-i \t-£ + ie t + Z-ie) (2P ■ q) 2 

| U(P 2 )f} U(P 1 )fS(t,0+mP2)^^j^)U(P 1 )g a (t,0 } • 

(51) 

Obviously the S a5fJjU , S a5 ^ ', e a5 ^ and e a5fJ,u expressions are obtained by the substitution IpTjl. 

At this stage, to square the amplitude, we need to calculate the following quantity, separately 
for the two charged processes 

T2 \rri |2 I rri rri* . rp rri* , I rri |2 /rr)\ 

— \J-DVNS\ + 1 DVNS 1 BH + 1 BH-L DVNS + I 1 BH\ , 

which simplifies in the neutral case, since it reduces \T DVNS \ 2 Eqs. (jl2J)-(jHJ and their axial 
counterparts are our final result and provide a description of the deeply virtual amplitude in 
the electroweak sector for charged and neutral currents. The result can be expressed in terms 
of a small set of parton distribution functions which can be easily related to generalized parton 
distributions, as in standard DVCS. 



8 Conclusions 

We have presented an application/extension of a method, which has been formulated in the past 
for the identification of the leading twist contributions to the parton amplitude in the generalized 
Bjorken region, to the electroweak case. We have considered the special case of a deeply virtual 
kinematics. We have focused our attention on processes initiated by neutrinos. From the the- 
oretical and experimental viewpoints the study of these processes is of interest, since very little 
is known of the neutrino interaction at intermediate energy in these more complex kinematical 
domains. 
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Abstract 

Neutrino scattering on nucleons in the regime of deeply virtual kinematics is studied both 
in the charged and the neutral electroweak sectors using a formalism developed by Bliimlein, 
Robaschik, Geyer and Collaborators for the analysis of the Virtual Compton amplitude in 
the generalized Bjorken region. We discuss the structure of the leading twist amplitudes of 
the process. 



Exclusive processes mediated by the weak force are an area of investigation which may gather a 
wide interest in the forthcoming years due to the various experimental proposals to detect neutrino 
oscillations at intermediate energy using neutrino factories and superbeams [1]. These proposals 
require a study of the neutrino-nucleon interaction over a wide range of energy starting from the 
elastic/quasi-elastic domain up to the deep inelastic scattering (DIS) region (see [2], [3], [4], [5], [6] 
for an overall overview). However, the discussion of the neutrino nucleon interaction has, so far, 
been confined either to the DIS region or to the form factor/nucleon resonance region, while the 
intermediate energy region, at this time, remains unexplored also theoretically. Clearly, to achieve 
a "continuos" description of the underlying strong interaction dynamics, from the resonant to the 
perturbative regime, will require considerable effort, since it is experimentally and theoretically 
difficult to disentangle a perturbative from a non-perturbative dynamics at intermediate energy, 
which appear to be superimposed. This is best exemplified - at least in the case of electromagnetic 
processes, such as Compton scattering - in the dependence of the intermediate energy description 
on the momentum transfer [7]. In this respect, the interaction of neutrinos with the constituents 
of the nucleon is no different, once the partonic structure of the target is resolved. From our 
viewpoint, the presence of such a gap in our knowledge well justifies any attempt to improve the 
current situation. 

Together with Amore, we have pointed out [8] that exclusive processes of DVCS-type (Deeply 
Virtual Compton Scattering) could be relevant also in the theoretical study of the exclusive neu- 
trino/nucleon interaction. Thanks to the presence of an on-shell photon emitted in the final state, 
this particle could be tagged together with the recoiling nucleon in a large underground detector 
in order to trigger on the process and exclude contamination from other backgrounds. With these 
motivations, a study of the vN — > uNj process has been performed in [8]. The process is mediated 
by a neutral current and is particularly clean since there is no Bethe-Heitler contribution. It has 
been termed Deeply Virtual Neutrino Scattering or DVNS and requires in its partonic description 
the electroweak analogue of the "non-forward parton distributions", previously introduced in the 
study of DVCS. 

In this work we extend that analysis and provide, in part, a generalization of those results to 
the charged current case. Our treatment, here, is purposely short. The method that we use for 
the study of the charged processes is based on the formalism of the non-local operator product 
expansion and the technique of the harmonic polynomials, which allows to classify the various 
contributions to the interaction in terms of operators of a definite geometrical twist [11]. We 
present here a classification of the leading twist amplitudes of the charged process while a detailed 
phenomenological analysis useful for future experimental searches will be given elsewhere. 

2 The Generalized Bjorken Region and DVCS 

Fig. 1 illustrates the process that we are going to study, where a neutrino of momentum / scatters 
off a nucleon of momentum Pi via a neutral or a charged current interaction; from the final state 
a photon and a nucleon emerge, of momenta q2 and P2 respectively, while the momenta of the 
final lepton is I'. We recall that Compton scattering has been investigated in the near past by 
several groups, since the original works [14, 15, 17]. A previous study of the Virtual Compton 
process in the generalized Bjorken region, of which DVCS is just a particular case, can be found 
in [16]. From the hadronic side, the description of the interaction proceeds via new constructs 
of the parton model termed generalized parton distributions (GPD) or also non-forward parton 
distributions. The kinematics for the study of GPD's is characterized by a deep virtuality of the 
exchanged photon in the initial interaction (v+p — > u+p + ^) ( Q 2 pa 2 GeV 2 ), with the final state 
photon kept on-shell; large energy of the hadronic system (W 2 > 6 GeV 2 ) above the resonance 
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Figure 1: Leading hand-bag diagrams for the process 



domain and small momentum transfers \t\ < 1 GeV 2 . In the electroweak case, photon emission 
can occur from the final state electron (in the case of charged current interactions) and provides 
an additional contribution to the virtual Compton amplitude. We choose symmetric defining 
momenta and use as independent variables the average of the hadron and gauge bosons momenta 

Pi,2 = PT^ ?i,2 = ?±y, (1) 
with A = P2 — Pi being the momentum transfer. Clearly 

A 2 

P ■ A = 0, P 2 = M 2 - — (2) 

and M is the nucleon mass. There are two scaling variables which are identified in the process, 
since 3 scalar products can grow large in the generalized Bjorken limit: g 2 , A • q, P ■ q. 

The momentum transfer A 2 is a small parameter in the process. Momentum asymmetries 
between the initial and the final state nucleon are measured by two scaling parameters, £ and rj, 
related to ratios of the former invariants 



^ = -2p-q ^=2p-q ® 

where £ is a variable of Bjorken type, expressed in terms of average momenta rather than 
nucleon and gauge bosons momenta. The standard Bjorken variable x = —q\j(2P\ ■ q\) is trivially 
related to £ in the t — limit and in the DVCS case r] = — £. 

Notice also that the parameter £ measures the ratio between the plus component of the mo- 
mentum transfer and the average momentum. 

£, therefore, parametrizes the large component of the momentum transfer A, which can be 
generically described as 

A = -2£P - A ± (4) 



where all the components of A^ are O (^y |A 2 

3 Bethe-Heitler Contributions 

Prior to embark on the discussion of the virtual Compton contribution, we quote the result for 
the Bethe-Heitler (BH) subprocess, which makes its first appearance in the charged current case, 
since a real photon can be radiated off the leg of the final state lepton. The amplitude of the BH 
contribution for a W + exchange is as follows 



T ^ = -\e\-^=-^u{l') 



[i — Ay + it 



2 



(5) 



where e is the polarization vector of the photon and 



T, 



w- 

BH 



9 9 _ m 
|e|- — — t=v{1) 



U(P2) 



r 



(Z- A) 2 + ie 



-7"(1 -7 5 ) 



:<(92) x 



(i^A 2 ) - i^(A 2 )) 7 a + (F 2 "(A 2 ) - F 2 d (A : 



' 2M 



^1) 



(6) 



for the W case, with D vS (qi)/(A 2 — + ie) being the propagator of the W's and F 12 the usual 
nucleon form factors (see also [8]). 

4 Structure of the Compton amplitude for charged and neu- 
tral currents 



Moving to the Compton amplitude for charged and neutral currents, this can be expressed in 
terms of the correlator of currents 



T^ = ij d A xe^(P 2 \T (jJ(:r/2) J™ ± > z °(-x/2)) \P ± ) , 
where for the charged and neutral currents we have the following expressions 



(7) 



j^(-x/2) = —L—^-x^-f^y + g z uA i^u(-x/2) + ^ d (-x/2)r(glv + g z dA i 5 )M-^m, 

A COS t>w 

J» w+ {-x/2) = ^f u (-x/2)Y(l - 7 5 RVM-*/2), 



r w {-x/2) 



2V2 



ii d {-x/2)^{l-^)U d M-x/2), 



Here we have chosen a simple representation of the flavour mixing matrix U* d = U ud = U du 
cos 9c, where 9c is the Cabibbo angle. 
The coefficients gy and g\ are 



(8) 



and 



7 L ' 2/i 

g u v = 2 + 3 sm °w 
1 2 

9dv = ~o + oSin 2 9 w 



gu g 5 g d 



gL 



g d A 



1 

"2 
1 

2 ' 



1 

3 



(9) 



(10) 



are the absolute values of the charges of the up and down quarks in units of the electron charge. 
A short computation gives 



(P 2 \T(j:(x/2)JZ°(-x/2))\P 1 ) = 
(P 2 $ u (x/2)g ulu S(x) ltM (gZ v + g z uAl ^ u {-x/2) 
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^ d {x/2) llx {g z dV + g z dA ^)S(-x)g dl ^ d ^/nPi), 

(11) 

{P 2 \T(j2(x/2)J^\-x/2))\P 1 ) = 
(P 2 \^(-x/2) lfi (l - ^)U ud S{-x) lu {-g d ) i/> d (x/2) + 
^ u {x/2) lv (g u ) S(x)^(l - ^)U ud M-*m\Pi), 

(12) 

(P 2 \T(rj(x/2)J^(-x/2))\P 1 } = 

(P 2 \ - ^ d {x/2)g dlv S{x) llx {\ - 1 b )U du ^ u {-x/2) + 

^ d (-x/2)^(l - 7 5 )5(-x)f/ du ^^/2)|Pi) , (13) 

where all the factors g/2\/2 and g/2cos9w, for semplicity, have been suppressed and we have 
defined 



S u (x) = S d (x) 



2tt 2 (x 2 — ie) 2 



(14) 



Using the following identities 



5 _ c P 5 _ • (i 



(15) 



we rewrite the correlators as 



Tf? = i / d A x 



- /IV 



2tt 2 (x 2 — ie) 



\9u9uV [ 



~9d9dv {S^ctvpO^ — ie^avpO^f^J + gd9dA 

(16) 



- /J.U 



rpW 
[IV 



T 

27r 2 (x 2 — ie) 2 



d 4 x 



iS, 



+ 



\Pl), 



iqx ajj 
■ i ;4 ° ^ du 

t ax- 



2n 2 (x 2 - ief 



~iS, av p (<5£ + O b £) - e, avP (Oi + Of)] \P X 



(17) 



(18) 



We have suppressed the x-dependence of the operators in the former equations. The relevant 
operators are denoted by 

Of (x/2, -x/2) = ^(x/2) 7 ^ (-x/2) + ^ a (-x/2)^ a (x/2), 
6f (x/2, -x/2) =j0 a (x/2) 7 5 7 ^a(-^/2) r f a (-x/2) 7 Vi(x/2), 
Of (x/2, -x/2) = ^(x/2) 7 > a (-x/2) - i/> a (-x/2)<fi/, a (x/2), 
Of '(x/2, -x/2) = V> a (*/2) 7 5 yVa(-z/2) + t(-x/2) 7 5 7^(x/2), 



Of d (x/2, -x/2) = g u ^ u (x/2)^ d (-x/2) + g^ u {-x/2)^ d (^l 
0%(x/2, -x/2) = ^ tt (x/2) 7 V^(-^/2) - ? A(^/2)7 5 7%(x/2), 



(19) 
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(20) 



and similar ones with u <-> d interchanged. 

We use isospin symmetry to relate flavour nondiagonal operators (Off) to flavour diagonal 
ones (Off) 



(p\O ud (x)\n) = (p\O ud (x)r~\n) = (p\ O ud (x),r~ \n) = (p\O uu (x)\p) - (p\O dd (x)\p) , 



\udi 



dd I 



(p\O ud (x)\n) = (n\O dd (x)\n) - (n\O uu (x)\n) , 
(n\O du (x)\p) = (p\O uu (x)\p) - (p\O dd (x)\p) , 



(n\O du (x)\p) = (n\O dd (x)\n) - (n\O uu (x)\n), 

where 

T ± = T X ± T V 

are isospin raising/lowering operators expressed in terms of Pauli matrices. 



(21) 
(22) 



5 Parameterization of nonforward matrix elements 

The extraction of the leading twist contribution to the handbag diagram is performed using the 
geometrical twist expansion, as developed in [10, 12, 13], adapted to our case. We set the twist-2 
expansions on the light cone (with x 2 = 0) and we choose the light-cone gauge to remove the 
gauge link 



(p^-kx^Mkxmy™- 2 = 

J Dze- ik ^F a ^( Zl , z 2 , P z ■ P 3 x\ P z ■ Pj)U(P 2 ) [Y - ikP/f } U(P l ) + 



J Dze- ik ^G a ^( Zl , z 2 , P t ■ P 3 x\ P % ■ P j )U(P 2 ) 
with < k < 1 a scalar parameter, with 

P z = PlZi+ P 2 z 2 , 

and 



M 



M 



U(Pi) , 
(23) 

(24) 



(P 2 \i; a (-kx)^rMkx)\Pi) tw - 2 = 

J Dze- ikix - Pz) F 5aM (z u z 2 , Pi ■ Pjx 2 , Pi ■ Pj)U(P 2 ) [ 7 V - ikP/-f 5 gt ] U(P X ) + 
J Dze- ik ^G 5a ^( Zl , z 2 , P t ■ P jX \ P t ■ Pj)U(P 2 h 5 



U(Pi). 
(25) 

The index (u) in the expressions of the distribution functions F, G has been introduced in order 
to distinguish them from the parameterization given in [10, 16], which are related to linear combi- 
nations of electromagnetic correlators. In the expressions above a is a flavour index and we have 



■5 



and P 2 , and the measure of integration is defined by [10] 



Dz = -dz x dz 2 9(1 - z ± ) 9(1 + z ± ) 9(1 - z 2 ) 9(1 + z 2 ). (26) 

In our parameterization of the correlators we are omitting the so called "trace-terms" (see 
ref. [16]), since these terms vanish on shell. In order to arrive at a partonic interpretation one 
introduces variables z + and Z- conjugated to 2P and A and defined as 

z + = l/2( Zl + z 2 ), 

Z_ = 1/2(32 -Zx), 

Dz = dz + dz_9(l + z + + z-)9(l + z + - *_)0(1 - z + + z_)0(l - z+ - z_). (27) 

In terms of these new variables P z = 2Pz + + Az_, which will be used below. 

At this stage, we can proceed to calculate the hadronic tensor by performing the x-space 
integrations. This will be illustrated in the case of the W + current, the others being similar. We 
define 



/ 
/ 
/ 
/ 



dx 



e iqx x a 



2tt 2 (x 2 - ie)- 



dx 4 



e %qx x a 



dx 4 



2-7T 2 (x 2 — ie)" 
2tt 2 (x 2 - ief 



■{P2\Sy,av()O a I Pi) — S* v , 
■(-P2I5 'pa,, pO ba \Pi) = Sf£, 



P2ViiavliO a ^\Pxl — 



dx 4 



2ix 2 (x 2 — ie) 



\P2\^^ av f30 bal3 \Pi) - e^l, 



(28) 



and introduce the variables 



(29) 



where (z) is now meant to denote both variables (z +1 zJ). The presence of a new variable Q 2 , 
compared to [10], is related to the fact that we are parameterizing each single bilinear covariant 
rather then linear combinations of them, as in the electromagnetic case. 
After some re-arrangements we get 



S;„ = g u j Dz 



I 



F a ^(z) 
(Ql + ie) 



(Ql + ie) 



G ai - V \z) 



{1- 



+ 



(Ql + ie) 



9u S Dz m 



+ ie) 



PzjjPzv PzpSlv PzvQ.il 9/-iv(Pz ' 0) 
+ Qvlix + %lv\ - [Pz^lv + Pzvlfi 

Pzfj,Pzu PzfiQu PzvQ.li 9fJ.v(Pz ' Q) 

+ 



+ 



-g^v M +Qv n/r + Qn 



M 



M 



+ 

"PA, 



M 



M 



6 



9d J Dz 



G a ^\z) 
(Qi + <ej 



-9,u U + 9,^— + %■ 



M 

M(Q 2 + ie) 



M 



^ M 



M 



(30) 



with an analogous expressions for that we omit, since it can be recovered by performing the 
substitutions 



(31) 



in (30). 

Similarly, for e a we get 



ft, J DzF'^(z- 1 



26 



(Ql + ie 



+ 



^ J DzG« v \z) 



(Q2 + k) W 



(Qi + ic] 

Q ^A g pP<f{ia xs g x A s ) 
M M(Q 2 + it) 

ta^A s pPq a (ia x5 q x A 5 y 
M M{Ql + ie) 



(32) 



The expression of e 5 ^ can be obtained in a similar way. 



To compute the tensor we need to include the following operators, which are related to 
the previous ones by g u , — > 1 



f ix > 
j dx > 

dx 4 

In this simple manipulation of (30) gives 

F a ^(z) 



2tt 2 (x 2 - 


ie) 






2tt 2 (x 2 - 


■ ^2 






2tt 2 (x 2 - 




e iqx x a 




2tt 2 (x 2 - 


ie) 



[P2\S 'nav pO al5 \Pi) = S° u , 

\P2\S ^ au pO bal3 \Pi) = s^i, 

7T2 ^ViictvpO"' |Pi) = , 



5a 



(33) 



v y 



(<5i + <c) 



/ 



F a ^(z) 



(Qi + ie 



{ [ Q\ivil + qul^ + q^lu] + \PzpT1v + Pzul^ 

PzpPzv Pz^qu PzvQn Qy,v\Pz ' q) 

Pz»lv + Pzvlp, 



(Qi + ie) 



{[-9^ + qui p. + q^lu] - 
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/ 



/ 



Dz 



Dz 



G a ^(z) 
(Qi + »e) 



G a( - V \z) 

(QI + 



9nu 



ia al3 Q n Afi ia^Af 



M 

M(Qj + ie) 



+ Qu- 



M 



M 



+ 



P i° vP Ap | P 



Zfl 



M 



M 



PzpPzv ~\~ PzfiQv ~\~ PzvQ.ii Q]xv\Pz ' Q) 



ia af3 q n A« ia^Ar- 



M 



+ q v - 



M 



ia uf3 A p 
M 



±7,11 „ , "T J 21; 



M 



ia a Pq a Ap 
~M(Q 2 2 + ie) 



PzfiPzv Pz^qv ~\~ Pzvq^i 9nv(Pz ' q) 



The expression of is obtained from (34) by the replacements (31). 
For the e a case, a re-arrangement of (32) gives 



§% = J DzF a ^{ 



(g2 + k) W 



g a 7 /3 - 



(Qi + ie) 



y DzG a( - v \z) 



+ 
+ 



(Qt + ie)\ 
a iaP s A s pPq a (ia x5 q x A 5 ) 



(g2 + k) W 



M M (Q 2 + ie) 

^(T^Ai P z V(icr A,5 g A A 5 ) 



M 



+ 



M(Q| + ie) 



+ 



Also in this case, the expression of the tensor is obtained by the replacements (31). 



M 



(34) 



(35) 



6 The partonic interpretation 

At a first sight, the functions F 5 ( u \ G 5 ^ do not have a simple partonic interpretation. 

To progress in this direction it is useful to perform the expansions of the propagators 



1 



1 



1 



Ql + ie 2(P ■ q) [z+ - £ + rjz- + ie] 
1 1 1 



Q\ + ie ' 2(P ■ q) [z+ + £ + rjz- - ie] 



(36) 



which are valid only in the asymptotic limit. In this limit only the large kinematical invariants 
and their (finite) ratios are kept. In this expansion the physical scaling variable £ appears quite 
naturally and one is led to introduce a new linear combination 



to obtain 



t = z + + r]z^ , 



1 1 



(37) 



1 



Q\ + ie 2(P -q)[t-S + ie]' 
1 1 1 



Q 2 2 + ie " 2(P-q)[t + Z-ie]' 
Analogously, we rewrite P z using the variables {t,Z-} 

P z = 2Pt + nz_ , 



(38) 



(39) 
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7T = A + 2£P . 



(40) 



This quantity is strictly related to Aj_, as given in (4). The dominant (large) components of 
the process are related to the collinear contributions, and in our calculation the contributions 
proportional to the vector n will be neglected. This, of course, introduces a violation of the 
transversality of the process of 0(Aj_/2P ■ q). 

Adopting the new variables {t, Z-} and the conjugate ones {P, 7r|, the relevant integrals that 
we need to "reduce" to a single (partonic) variable are contained in the expressions 
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2P-qJ 



Dz 



(Ql + ie) 



H(t + £z^zj 
(t-£ + ie) 



2PH + t^z. 



(Ql + ie) 

H(t + £z-,z-) 
{t-t + ief 

+7iVz 2 _ + (q^ + q v ^) z_ + (2PV + 2PV) tz. 



Dz- 



4P»P u t 2 + (2q^P v + 2<f P") t - g^(q ■ P z ) 



(41) 



Here H(z + , zJ) is a generic symbol for any of the functions. We have similar expressions for the 
integrals depending on the momenta Q2. 

The integration over the z_ variable in the integrals shown above is performed by introducing 
a suitable spectral representation of the function H(t, +£z~, Z-). As shown in [10], we can classify 
these representations by the n — 0, 1, powers of the variable z_ , 



h n (t/r,0= j dz.z. n h(- + iz.,zJ). 
With the help of this relation one obtains 



(42) 



If 1 f 1 dr 

H n (t,0 = -J dz^ n H(t + ^,z.) = -J -r n h n (t/r,0 



(43) 



The result of this manipulation is to generate single- valued distribution amplitudes from double- 
valued ones. In the single-valued distributions h n (t,£) the new scaling variables t and £ have a 
partonic interpretation. £ measures the asymmetry between the momenta of the initial and final 
states, while it can be checked that the support of the variable t is the interval [—1,1]. The twist-2 
part of the Compton amplitude is related only to the n = moment of Z-. Before performing the 
z_ integration in each integral of Eq. (41) using Eq. (43) - a typical example is Hq^) - we reduce 
such integrals to the sum of two terms using the identity 



dt- 



t" 



dt- 



t 



m—l 



H n (t,0--h n (t,0 



(44) 



-1 (t±£Tie) 2 " v ' s/ J-i (t±£Tie) 
As shown in [16], after the Z- integration, the integrals in (41) can be re- written in the form 

1 



H Ql (0 



2P ■ q J-i (t - £ + ie) ' 
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mo 



(2P.q) 



— I 1 

■ ■ q) 2 J-i 



dt 



2H (t,0-h (t,0 



(2P ■ q) 2 J-i 



dt 



(t - £ + ie) 
H (t,O-h (t,O] 



Up^P u t 



{(2<f + 2q v P> x - g^2q ■ P)} + 0(nV), 



(45) 



where, again, we are neglecting contributions from the terms proportional to tt^, subleading in 
the deeply virtual limit. The quantities that actually have a strict partonic interpretation are the 
ho(t,C) functions, as argued in ref. [19]. The identification of the leading twist contributions is 
performed exactly as in [10]. We use a suitable form of the polarization vectors (for the gauge 
bosons) to generate the helicity components of the amplitudes and perform the asymptotic (DVCS) 



limit in order to identify the leading terms. Terms of 0(1/ \J2P ■ q) are suppressed and are not 
kept into account. Below we will show only the tensor structures which survive after this limit. 



7 Organizing the Compton amplitudes 



In order to give a more compact expression for the amplitudes of our processes we define 

g T,u = _ g ,u + 

a(t) 



Jq 7 P) + OTP)' 
9u 9d 

(t-S + ie) (t + £- ie) ' 

9u 9d 

(t-S + ie) (t + £- ie) • 



(46) 



Calculating all the integrals in the Eqs. (30), (34) and (32), (35), we rewrite the expressions of 
the amplitudes as follows 



T% + = iU ud U(P 2 ) [i (S; u + S%) + e% + e% - i (fij, + S%) - e%, - e£] U(P l ) , 



^% — —iU du U(P 2 ) i (S™ v + S^fj + £^ v + e^ v — i (s^, + S^fj — — e 
T* = iU(P 2 ) [g u g uV (S; u - ie£) - g u g uA (§% - i~e%) 

-g d g d v (sj, - ie%) + g d9dA (S% - ie%)] U(P 1 ) , 
where, suppressing all the subleading terms in the tensor structures, we get 



U(Pi) , 



(47) 



U(P 2 )S a ^U(P l 



1 gT^V 

, dta{t) 2p—q 



u{p 2 )i u^fsit, o + u(p 2 ) (i^j^mp.mt, o 



+ 



J 1 dtp(t) 



pti pu 

(p-qj< 



U{P 2 )i UiPjtfSfrt) + U(P 2 )(i^j^)U(P 1 )tg a (t,0 



(48) 



while for the e afiu expression we obtain 



U(P 2 )e^U(P 1 )=e^0^ j^dtp(t) [u(P 2 )i U(P 1 )fS(t,0 



+ 



U(P2)(i^^)U(P 1 )g a (t,0 



(49) 
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U(P 2 )S a ^U(P 1 



dt 



t — £ + it t + £ — it 



x 



9 



[2P-q 

/>( 



u{P2)i rn' : )/;;u.a ■ u(p 2 )(i a 2 



M 



mp.mt^) 



+ 



i 



t-t + ie t + t-iej (P ■ qf 



u{PM u(PM(t,o + u(p 2 ){i^^-)u(p 1 )mu) 



(50) 



U{P 2 )e a ^U(P 1 ] 



dt 



t — £ + it t + £ — it 



(2P • qf 



U(P 2 )i U(P 1 )fS(t : + U(P 2 )(i 



M 



)U(P 1 )g a (t,0 



(51) 



Obviously the S 115 ^, S a5fMU , e a5 ^ u and e a5 ^ expressions are obtained by the substitution (31). 

At this stage, to square the amplitude, we need to calculate the following quantity, separately 
for the two charged processes 



T2 \rri 1 2 I rri rTi* , rji rri* , I rri |2 

— \J-DVNS\ + -LDVNS-LbH + J-BH-I-DVNS + \ I BH\ , 



(52) 



which simplifies in the neutral case, since it reduces \T DVNS \ 2 [8]. Eqs. (47)-(51) and their axial 
counterparts are our final result and provide a description of the deeply virtual amplitude in 
the electroweak sector for charged and neutral currents. The result can be expressed in terms 
of a small set of parton distribution functions which can be easily related to generalized parton 
distributions, as in standard DVCS. 



8 Conclusions 



We have presented an application/extension of a method, which has been formulated in the past 
for the identification of the leading twist contributions to the parton amplitude in the generalized 
Bjorken region, to the electroweak case. We have considered the special case of a deeply virtual 
kinematics. We have focused our attention on processes initiated by neutrinos. From the the- 
oretical and experimental viewpoints the study of these processes is of interest, since very little 
is known of the neutrino interaction at intermediate energy in these more complex kinematical 
domains. 
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